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Abstract. A well-known potential dark matter sig- 
nature is emission of GeV - TeV neutrinos from 
annihilation of neutralinos gravitationally bound to 
massive objects. We present results from recent 
searches for high energy neutrino emission from 
the Sun with AMANDA, in all cases revealing no 
significant excess. We show limits on both neutralino- 
induced muon flux from the Sun and neutralino- 
nucleon cross section, comparing them with recent 
IceCube results. Particularly, our limits on spin- 
dependent cross section are much better than those 
obtained in direct detection experiments, allowing 
AMANDA and other neutrino telescopes to search a 
complementary portion of MSSM parameter space. 
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I. Introduction 

Weakly interacting massive particles (WIMPs) with 
electroweak scale masses are currently a favored expla- 
nation of the missing mass in the universe. Such particles 
must either be stable or have a lifetime comparable to 
the age of the universe, and they would interact with 
baryonic matter gravitationally and through weak inter- 
actions. The minimal supersymmetric standard model 
(MSSM) provides a natural candidate, the lightest neu- 
tralino A large range of potential neutralino masses 
exists, with a lower bound on the mass of the lightest 
neutralino of 47 GeV imposed by accelerator-based 
analyses [2], while predictions based on the inferred dark 
matter density suggest masses up to several TeV 0. 

Searches for neutralino dark matter include direct 
searches for nuclear recoils from weak interaction of 
neutralinos with matter J4), and indirect searches 
for standard model particles produced by neutralino 
annihilation. Particularly, a fraction of neutralinos inter- 
acting with massive objects would become gravitation- 
ally bound and accumulate in the center. If neutralinos 
comprise dark matter, enough should accumulate and an- 
nihilate to produce an observable neutrino flux. Searches 
for a high energy neutrino beam from the center of 
the Earth and the Sun 0, 0, 0, (H, flZQ have 
yielded negative results. Observations of a cosmic ray 
electron-positron excess by ATIC fl21, PPB-BETS fT51 . 
Fermi [14|, and HESS 1151 . along with the anomalous 
cosmic ray positron fraction reported by PAMELA fl6l . 
could be interpreted as an indirect signal of dark matter 
annihilation in our galaxy ifTTl . 



Here we present searches for a flux of GeV-TeV 
neutrinos from the Sun using AMANDA. We improve on 
the sensitivity of the previous AMANDA analysis IfTTl 
significantly and extend the latest results from IceCube 
[7 1 to lower neutralino masses. We observe no neutralino 
annihilation signal and report limits on the neutrino- 
induced muon flux from the Sun and the resulting limits 
on neutralino-proton spin-dependent cross section. 

II. Neutrino Detection with AMANDA 

The detection of neutrino fluxes above ~ 50 GeV is a 
major goal of the Antarctic Muon And Neutrino Detector 
Array (AMANDA). AMANDA consists of 677 optical 
modules embedded 1500 m to 2000 m deep in the ice 
sheet at the South Pole, arranged in 19 vertical strings 
and occupying a volume of ~ 0.02 km 3 . Each module 
contains a 20 cm diameter photomultiplier tube (PMT) 
optically coupled to an outer glass pressure sphere. 
PMT pulses ("hits") from incident Cherenkov light are 
propagated to surface electronics and are recorded as 
an event when 6-7 hits on any one string or 24 total 
hits occur within 2.5 /is. The vast majority of the 
O(10 9 ) events recorded each year are downgoing muons 
produced by cosmic ray air showers in the atmosphere 
above the South Pole. Relativistic charged leptons pro- 
duced near the detector via charged-current neutrino 
interactions similarly trigger the detector, with several 
thousand atmospheric neutrino induced muon events 
recorded per year. The hit leading edge times, along 
with the known AMANDA geometry and ice properties 
lfl8l . allow reconstruction of muon tracks with median 
accuracy 1.5° - 2.5°, dependent on zenith angle. 

AMANDA operated in standalone from 2000-2006 
and is currently a subdetector of the much larger (~ 
km 3 ) IceCube Neutrino Observatory |fL9ll , scheduled 
for completion in 2011. The optical module density of 
AMANDA is much higher than that of IceCube, making 
AMANDA more efficient for low-energy muons (< 300 
GeV) which emit less Cherenkov light. 

III. Data Selection and Methods 

We describe two separate searches for Solar neu- 
tralinos in this proceeding. First, we present a search 
using a large data sample from 2000-2006 prepared 
for a high energy extraterrestrial point source search 
l20l . l2ll . We also present a search using data from 
2001-2003, optimized to retain low energy events l22l . 
Both analyses are done in two stages; first, neutrino 
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induced muon events are isolated from the much larger 
background of downgoing muons, then a search method 
is used to test for an excess at the location of the Sun. 

A. Data Selection 

While the Sun is above the horizon, neutrino-induced 
muons from the Sun are masked by the much larger 
background of downgoing cosmic ray muons; thus, we 
select data during the period when the Sun is below the 
horizon (Mar. 21 - Sept. 21), resulting in 953 days live- 
time from 2000-2006 and 384 days from 2001-2003. In 
both analyses, neutrino events are isolated by selecting 
well reconstructed upgoing muon tracks. Events are first 
reconstructed with fast pattern matching algorithms, and 
events with zenith angles 9 < 80° (9 < 70° for the 
2001-2003 analysis) are discarded, eliminating the vast 
majority of downgoing muons. The remaining events 
are reconstructed with a more computationally intensive 
maximum-likelihood reconstruction |23| accurate to 1.5° 
- 2.5°, and again events with 9 < 80° are discarded. 

O(10 6 ) misreconstructed downgoing muon events re- 
main per year, and these are reduced by cuts on 
track quality parameters such as track angular uncer- 
tainty [24 1, the smoothness (evenness) of hits along 
the track l23l . and the likelihood difference between 
the maximum-likelihood track and a forced downgoing 
likelihood fit using the zenith distribution of downgoing 
muons as a prior [23 1. For the 2000-2006 analysis, 6595 
events remain after quality cuts, dominantly atmospheric 
neutrinos ll20l . reduced to 4665 events by requiring dates 
when the Sun is below the horizon. Zenith distributions 
from 2000-2006 are shown in figure Q] 

We consider neutralino masses from 100 GeV to 5 
TeV and two extreme annihilation channels: W + W~ 
(t + t~ for 50 GeV) and bb, which produce high and 
low energy neutrino spectra, respectively, relative to the 
neutralino mass. The fraction of signal events retained 
depends on the neutrino energy spectrum and varies from 
17% for a 5 TeV neutralino mass and W + W~ channel 
to 1% for 100 GeV and bb channel, relative to trigger 
level, in the 2000-2006 analysis. 

The 2001-2003 analysis is a dedicated neutralino 
search, unlike the 2000-2006 analysis, and more con- 
sideration is given to low energy events. Twelve event 
observables are considered, and selection criteria based 
on these observables are optimized separately for three 
signal classes, dependent on neutralino mass and anni- 
hilation channel, to maximize retention of signal events. 
The signal classes are shown below along with the 
number of events passing selection criteria. 



Class Channel 



Final Events 



W + W~ 
bb 

w + w~ 
bb 

T + T - 
bb 



5 TeV 
5 TeV 



250 GeV 
500 GeV 

100 GeV 
250 GeV 

50 GeV 
50, 100 GeV 



670 



504 



398 



10' 



10" 



10' r 



10 = 



£ 

w 10 4 



IO'e 



10 2 E 



10 



Downgoing Muons 



Misreconstructed Muons 



Filtered Data 



Atmospheric Neutrinos I 



*'(l-i'rt-iy«, 1Jpli i 

Final Sample i* 



Neutralino Signal 



1 i! 



-0.5 





cose 



0.5 



Fig. 1 . Reconstructed zenith angles of data (circles) at trigger level, 
filter level (8 < 80°), and final cut level, true (fine dotted) and 
reconstructed (solid) zenith angles of CORSIKA | 31| downgoing muon 
simulation at trigger level, reconstructed zenith angles of ANIS [28] 
atmospheric neutrino simulation at trigger level, filter level, and final 
cut level (dashed), and reconstructed zenith angles of a neutrino signal 
from the Sun at final cut level (dash-dotted). 



The selection is more efficient than the 2000-2006 
analysis, with 21% of signal retained for 5 TeV, W + W~ 
channel, to 4% for 100 GeV, bb channel. The 2001- 
2003 analysis additionally considers 50 GeV neutralino 
masses, with a signal efficiency of l%-3%. 

B. Search Method 

Both analyses use maximum-likelihood methods 
to search for an excess of events near the location of 
the Sun. The data is modeled as a mixture of n s signal 
events from the Sun and background events from both 
atmospheric neutrinos and misreconstructed downgoing 
muons. The signal likelihood for the i th event is 



Si = 



27TCT 2 



(1) 



where is the space angle difference between the event 
and the Sun, and <7, is the event angular uncertainty 
11241 . The background likelihood Bi is obtained from 
the zenith distribution of off-source data. The full-data 
likelihood over all N data events is 



N 

C = P{Data\n s ) = J] + (1 - ^)B t ) 



(2) 
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and is numerically maximized to find the best fit event 



excess n s . The likelihood ratio —2 log 



C(0) 
C{fi s 



is approxi- 



mately x 2 distributed and provides a measure of signif- 
icance. Event upper limits are set from this likelihood 
using the Feldman-Cousins unified construction [26|. 

C. Signal Simulation and Systematic Uncertainties 

Neutrino energy distributions at Earth from neutralino 
annihilation in the Sun are generated by DarkSUSY 
ll27l . For the 2000-2006 analysis, neutrino events are 
generated with ANIS |28|, with muons propagated using 
MMC [29 1, then reweighted to the energy distributions 
described above. For 2001-2003, the DarkSUSY energy 
distributions are sampled by WimpSimp |30|, and muons 
are propagated with MMC. 

Uncertainties in our signal simulation are dominated 
by uncertainties in optical module sensitivity and photon 
propagation in ice. These uncertainties are constrained 
by comparing the trigger rate of CORSIKA [31] down- 
going muon simulation using various hadronic models 
with the observed AMANDA trigger rate. The effect on 
signal prediction is measured by shifting the simulated 
optical module efficiency by these constraints and is 
10% for m x = 5 TeV, W+W~ channel, to 21% for 
m x = 100 GeV, bb channel. Other sources of uncertainty 
include event selection (4%-8%) and uncertainty in neu- 
trino mixing angles (5%). For the 2000-2006 analysis, 
uncertainties total 13%-24% and are included in the 
limit calculation using the method of Conrad et al. lf32l 
as modified by Hill 03 • Uncertainties for 2001-2003 
total 23%-38% and are included in the limits assuming 
the worst case. 

IV. Results 

The search methods are applied to the final data, and 
both analyses reveal no significant excess of neutrino- 
induced muons from the direction of the Sun. A Sun- 
centered significance skymap from the 2000-2006 anal- 
ysis (figure^ shows a 0.8er deficit from the direction of 
the Sun. For the 2001-2003 analysis, a deficit of events 
is observed in classes A and C, and a small excess is 
seen in class B. Each excess or deficit is within the la 
range of background fluctuations. 

Upper limits on the neutralino annihilation rate in the 
Sun are calculated from the event upper limit (iqq by 



N AP T L V eff 



dN„ 



-dE 



(3) 



,o dE 

where R is the Earth-Sun radius, Na is the Avogadro 
constant, p is the density of the detector medium, Tl 
is the livetime, and ovjv is the neutrino-nucleon cross 
section. The muon neutrino energy spectrum ^g- for a 
given annihilation channel is obtained from DarkSUSY 
and includes absorption and oscillation effects from tran- 
sit through the Sun and to Earth. The energy-averaged 
effective volume V e ff is obtained from simulation. Lim- 
its on muon flux are given by 



and limits on neutralino-proton cross section are calcu- 
lated according to IT341 . These quantities are tabulated in 
tableUfor the more restrictive of the two analyses. Muon 
flux limits, assuming a 1 GeV threshold on muon energy, 
and spin-dependent cross section limits are shown in 
figure [3] for both analyses. 

V. Discussion 

These limits extend the latest IceCube limits to lower 
neutralino masses and are now beginning to exclude 
neutralino spin-dependent cross sections allowed by 
direct detection experiments (figure [3}. A 1000-fold 
improvement over current direct-detection limits J4), 
does not significantly constrain allowed spin-dependent 
cross sections; thus, neutrino telescopes will continue to 
observe a complementary portion of MSSM parameter 
space over the next several years. IceCube is currently 
operating with 59 strings and will contain 86 strings 
when complete in 2011. The DeepCore extension to 
IceCube IT351 . six strings with tighter string spacing (72 
m), tighter optical module spacing (7 m), and higher 
PMT quantum efficiency, will be complete in 2010. 
DeepCore will significantly enhance the sensitivity of 
IceCube to low energy muons, extending the reach of 
IceCube to lower neutralino masses. 
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Fig. 2. Sun-centered skymap of event excesses from the 2000-2006 analysis. 
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TABLE I 

Effective volume, event upper limit, and preliminary limits on neutralino annihilation rate in the Sun, 
neutrino-induced muon flux from the sun, and spin-independent and spin-dependent neutralino-proton cross section 

for a range of neutralino masses, including systematics. 




Fig. 3. Preliminary limits on neutrino-induced muon flux from the Sun (left) along with limits from IceCube |7 |, BAKSAN |8], MACRO (5), 
and Super-K 1 10], and limits on spin-dependent neutralino-proton cross section (right) along with limits from CDMS |4|, IceCube |7 |, Super-K 
flOl , KIMS (36), and COUPP (37) • Th e 

green shaded area represents models from a scan of MSSM parameter space not excluded by the 
spin-independent cross section limits of CDMS [4| and XENON |5 |, and the blue shaded area represents allowed models if spin-independent 
limits are tightened by a factor of 1000. The projected sensitivity of 10 years operation of IceCube with DeepCore is shown in both figures. 



